function Rx = Rotx(theta)
    Rx = zeros(3,3);
    Rx = [1               0      0;
          0        cos(theta) -sin(theta);
          0        sin(theta)  cos(theta)];
end